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Abstract. We present Infrared Space Observatory (ISO) 
spectra of fourteen isolated Herbig Ae/Be (HAEBE) stars, 
to study the characteristics of their circumstellar dust. 
These spectra show large star-to-star differences, in the 
emission features of both carbon-rich and oxygen- rich dust 
grains. The IR spectra were combined with photometric 
data ranging from the UV through the optical into the 
sub-mm region. We defined two key groups, based upon 
the spectral shape of the infrared region. The derived re- 
sults can be summarized as follows: (1) the continuum of 
the IR to sub-mm region of all stars can be reconstructed 
by the sum of a power-law and a cool component, which 
can be represented by a black body. Possible locations 
for these components are an optically thick, geometrically 
thin disc (power-law component) and an optically thin 
flared region (black body); (2) all stars have a substan- 
tial amount of cold dust around them, independent of the 
amount of mid-IR excess they show; (3) also the near-IR 
excess is unrelated to the mid-IR excess, indicating dif- 
ferent composition/location of the emitting material; (4) 
remarkably, some sources lack the silicate bands; (5) apart 
from amorphous silicates, we find evidence for crystalline 
silicates in several stars, some of which are new detections; 
(6) PAH bands are present in at least 50% of our sample, 
and their appearance is slightly different from PAHs in the 
ISM; (7) PAH bands are, with one exception, not present 
in sources which only show a power-law continuum in the 
IR; their presence is unrelated to the presence of the sil- 
icate bands; (8) the dust in HAEBE stars shows strong 
evidence for coagulation; this dust processing is unrelated 
to any of the central star properties (such as age, spectral 
type and activity). 
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1. Introduction 

A circumstellar (CS) disc is expected to be a natural 
byproduct of the star forming process (e.g. Shu et al. 
1987). This theoretical expectation has obtained wide sup- 
port from optical (e.g. McCaughrean & O'Dell 1996), in- 
frared (e.g. Marsh et al. 1995) and millimetre observations 
of young stars (e.g. Mannings & Sargent 1997). The CS 
disc is expected and observed to gradually disappear, but 
remnants are still found around several Main-Sequence 
(MS) stars, such as Vega (Aumann et al. 1984). Earlier 
and recent modelling of T Tauri discs have shown that the 
most successful models are flaring passive discs (Kenyon 
& Hartmann 1987, Chiang & Goldreich 1997). 

Herbig Ae/Be stars (hereafter HAEBEs), first de- 
scribed as a group by Herbig (1960), are believed to be 
the more massive analogues of T Tauri stars. They are 
seen as the progenitors of Vega-type stars (for recent re- 
views, see Waters & Waelkens 1998; Natta et al. 1998). 
They are characterized by large IR excesses due to ther- 
mal re-emission of CS dust, show emission lines in their 
spectrum due to CS gas and have masses between 2 and 8 
M© (Herbig 1960). Infrared spectroscopy offers a unique 
opportunity to scrutinize the composition and character- 
istics of their CS dust. Recent ISO (Kessler et al. 1996) 
studies have revealed a large variety in the properties of 
the dust around HAEBEs, from which it became clear 
that their dust is significantly different from that in the 
interstellar medium (Waelkens et al. 1996, Malfait et al. 
1998a, Malfait et al. 1999b, van den Ancker et al. 1999). 



This paper is one in a series of papers based upon 
ISO-SWS observations of HAEBE stars. In this study, we 
compiled a set of data which include, next to the ISO 
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Fig. 1. ISO spectra of the 14 sample stars, superimposed on their spectral energy distributions. Crosses: observations; 
full line through the optical data: Kurucz model atmosphere; other full line: ISO-SWS/LWS observations; arrows 
indicate upper limits. The data are normalized to the V band. 



spectra, also UV, optical, IR and sub-mm photometry of 
a large sample of isolated HAEBE stars. A similar study 
was already presented by Sylvester et al. (1996) for a 
sample of Vega-like systems. Their ground-based observa- 
tions in the IR with UKIRT are restricted to 2 ranges: 
7.5-13.5 /im and 15.8-23.9 ^m. Some of their sources 



(HD135344, HD139614, HD142666, HD144432, HD169142 
and 51 Oph) are also part of our HAEBE sample, and it is 
interesting to compare their results with ours. In this pa- 
per we give an overview of the IR features in our sample, 
together with a description of the Spectral Energy Distri- 
butions (SED) and we propose a global model to explain 
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Fig. 2. The ISO-SWS spectra of our programme stars. Together with HD100546, we also show the spectrum of comet 
Halc-Bopp (Crovisier et al. 1997) for comparison. 



the SEDs. In Sect. 2, we describe our sample stars and 
their observations. We also present the SEDs (see Fig. |l]) 
and indicate observational trends. ISO-SWS spectra and 
an inventory of solid state and PAH bands are shown in 
Sect. 3, where the individual sources are discussed as well. 
In Sect. 4 we propose a global model, and discuss grain 
processing. Our conclusions are summarized in Sect. 5. In 
a forthcoming paper, detailed radiative transfer models of 



some of the sources will be presented (Bouwman, Meeus 
& Waters, in preparation). 

2. Targets and Observations 

The objects we selected are fourteen so-called isolated 
HAEBE stars. These HAEBE stars are not located inside 
a star-forming region, but show all the other characteris- 
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Table 1. Parameters of the fourteen programme stars. The groups are defined below (see Sect. 2.1). 







(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


Group 


Object 


Spectral 


T,,ff 


log g 


d 


log (Age) 


E[B - V]cs 


F800(imd 


Aonset 


Lir/L* 






type 


(KJ 




(P c ) 


\J r ) 




(mjy pc 2 ) 


(/mi) 






AB Aur 


B9/A0Ve 


9750 


5.0 


144 


6.3 


0.07 


11 10 6 


1.1 


0.48 


la 


HD 100546 


B9Ve 


11000 


4.5 


103 


> 7.0 


0.00 


12 10 6 


1.2 


0.51 




HD 142527 


F7IIIe 


6250 


4.0 


200 


5.0 


0.08 


166 10 6 


1.1 


1.06 




HD179218 


B9e 


10000 


5.0 


240 


5.0 


0.06 


11 10 6 


1.7 


0.62 




HD100453 


A9Ve 


7500 


4.5 






0.06 




1.2 


0.54 


lb 


HD 135344 


F4Ve 


6750 


4.5 


84 




0.00 


4 10 6 


1.1 


0.44 




HD139614 


A7Ve 


8000 


4.5 


151 




0.01 


14 10 6 


1.5 


0.39 




HD169142 


A5Ve 


10500 


4.5 


145 




0.00 


12 10 6 


1.6 


0.10 




HD104237 


A4Ve 


10500 


4.5 


116 


6.3 


0.25 


3 10 6 


1.2 


0.13 




HD 142666 


A8Ve 


8500 


4.5 


116 




0.40 


4 10 6 


1.1 


0.28 


Ila 


HD 144432 


A9Ve 


8000 


4.5 


> 200 




0.05 


> 4 10 6 


1.1 


0.26 




HD150193 


AlVe 


10000 


4.0 


150 


> 6.3 


0.30 


3 10 6 


1.2 


0.15 




HD163296 


A3Ve 


10500 


4.0 


122 


6.6 


0.02 


34 10 6 


1.2 


0.16 




51 Oph 


AOVe 


10000 


4.0 


131 


5.5 


0.03 


<.8 10 6 


2.3 


< 0.024 



References: (1): Malfait et al. (1998b), Dunkin et al. (1997), Gray & Corbally (1998); (2), (3) and (6): Malfait et al. (1998b); 
(4) and (5): van den Ancker et al. (1999); (7) and (9): this study, based upon sub-mm measurements by Sylvester et al. (1996), 
Mannings & Sargent (1997), Walker & Butner (1995), Henning et al. (1998) and Mannings & Sargent (2000); (8): this study. 



tics of a HAEBE star and are presumably the somewhat 
more evolved members of the HAEBE group. These ob- 
jects are best suited for our purpose, which is to discuss the 
evolution of the CS disc, and offer the additional advan- 
tage that the spectra are not strongly affected by loosely 
bound remnants of the star formation process, but only 
show the emission of the disc. In Table g we list the pro- 
gramme stars and among others their main parameters: 
spectral type, effective temperature and log g. 

The sources have been observed with the ISO Short 
Wavelength Spectrometer (SWS; de Graauw et al. 1996) 
in mode AOT1. The spectra cover an interval from 2 to 
45 /zm. Some stars have also been observed with ISO- 
LWS (Clegg et al. 1996), which covers a range from 45 
to 200 /j,m. These spectra were discussed by Malfait et al. 
(1998a, 1999a, b). In this study we will concentrate on the 
SWS data. The spectra were reduced in a standard way us- 
ing the ISO-SWS Interactive Analysis (IA) tool containing 
pipeline processing steps of OLP version 8.5, and the ISO 
Spectral Analysis Package (ISAP version 1.6a). In Fig. | 
the reduced SWS spectra are shown. For some sources 
(HD135344, HD150193 and 51 Oph), the ratio signal-to- 
noise is so low at longer wavelengths that we had to leave 
out the part longwards of 28 /im. 

2.1. Spectral Energy Distributions 

We also collected photometric data in the literature (Mal- 
fait et al. 1998b and references therein; sub-mm data from 
Sylvester et al. 1996, Mannings & Sargent 1997, 2000; 



Walker & Butner 1995, and Henning et al. 1998), and com- 
posed for each star an SED, ranging from the UV until the 
sub-mm region. In Fig. [l], we show for each of the fourteen 
sources its SED, combined with their respective ISO spec- 
trum. An appropriate Kurucz (1993) model atmosphere 
was fitted through the optical data, representing the pho- 
tospherical contribution; it emphasizes the shape and the 
amount of the excess in the IR and sub-mm region. The 
shortest wavelength at which an excess is discernible is 
listed in Table [l] as A on sct, with an uncertainty of 0.2 /im. 
Also shown in Table [l] is the derived fractional luminosity 
of the dust, Lir,/L*, which is the ratio of the energy radi- 
ated by the dust to the stellar luminosity. This ratio was 
calculated as follows: first, we converted the data into the 
Fa versus A scale. Then we integrated the Kurucz model 
over its entire wavelength range to calculate L*. To ob- 
tain Lir, we first subtracted the Kurucz model from the 
observations, and then integrated this curve longwards of 
Aonset- Sylvester et al. (1996) have already calculated this 
ratio for six of our sample stars, and their results agree 
very well with four of our stars, while they agree less well 
for HD135344 (0.64) and HD1444432 (0.48). The values 
we obtained for Lir/L* are consistent with a passive re- 
processing disc, except for the star HD142527 (Lir/L* = 
1.06). 

The dust continuum behaves very differently from 
source to source, especially in the mid-IR (15 to 45 /mi). 
In some stars it is rising, in other stars it is rather flat or 
even descending. Also the strength of the dust continuum 
in these objects is very diverse: the 12 /tm excess ranges 
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between 3.5 and 7 magnitudes, the 60 /Ltm excess ranges 
between 4.5 and 12 magnitudes and the 1.3 mm excess 
between 10 and 13 magnitudes. A second important ob- 
servational fact is the strong variation of the strength of 
the 10 /im silicate feature from one object to another. In 
some objects (e.g. HD144432) this feature is very strong, 
in others (e.g. 51 Oph) it is less so, and in several objects 
(e.g. HD169142) it is even absent. 

Notwithstanding these large differences, the overall 
structure of the dust discs seems to be similar. It is pos- 
sible to decompose the spectra into at maximum three 
components: a power-law, a black body (BB) and the 
solid state bands. As a first step, we fitted the IR con- 
tinuum of the stars showing a flat continuum with a 
power-law. Actually, the determination of the continuum 
is non-trivial and should be taken with some caution. Af- 
ter some experiments, we found that the continuum can 
be best determined by plotting the spectrum as log F\ 
versus log A. As an example, we show in Fig. ||, up- 
per panel, how the continuum determination was done 
for HD150193. It is surprising to see that for at least 
six sources (HD104237, HD142666, HD144432, HD150193, 
HD 163296 and 51 Oph), the continuum can be fitted very 
well with a power-law. We have to remark here that we did 
not remove the photospheric component, since it is only a 
negligible (< 10 %) fraction of the total flux. Only for one 
source, 51 Oph, the dust component is less dominant, and 
the photospheric contribution to the total flux in the IR 
is much more important. We therefore did not determine 
the power-law continuum for this source. 

We then proceeded to fit the sources with a rising mid- 
IR continuum, assuming we could apply a similar power- 
law continuum fit for these sources. This assumption is 
supported by the similarity in the near-IR region for all 
our sample stars. For these sources, such as HD179218 
(see Fig. |[ lower panel), an additional black bodyQ (BB) 
on top of a power-law is needed to fit the continuum. Re- 
markably, with just these two components (a power-law 
and a BB), the continuum of all sources can be fitted. The 
slope of the power-law fits are listed in Table ||, together 
with the average temperature of the respective BBs, when 
needed. In the sub-mm region, we observe a turn-down in 
the slope of the continuum. This is because sub-mm wave- 
lengths are longwards of the peak of the black body for 
the coolest grains in the disc (thus we observe a Rayleigh- 
Jeans tail in the SED). 

We have accordingly classified the sample: group I 
contains sources for which the continuum can be recon- 
structed by a power-law and a black body, and group 
II sources only need a power-law to fit their continuum. 
The groups can be further subdivided according to the 

1 This BB component is usually a single-temperature BB, 
but it also can be a set of typically 2-3 BBs with different 
temperatures, to which we refer to as 'the' BB component for 
convenience, as an analogue to the power-law continuum. 



HD150193 




too 
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Fig. 3. The determination of the continuum for the 
sources HD150193 (group Ha) and HD179218 (group la) 
in log Fa versus log A space. Full line: SWS spectrum, 
straight dotted line: power-law continuum fit; curved dot- 
ted line: sum of a power-law and a BB (T ~ 190 K) con- 
tinuum; asterisks: IRAS colour-corrected fluxes. 



presence or absence of solid state bands. In group Ia/IIa 
sources solid state bands are present, while group Ib/IIb 
sources are without solid state bands. Our sample does not 
contain a source which would fit in group lib, i.e. there is 
no star in our sample with a pure power-law continuum 
which lacks solid state bands. This can be an observational 
selection effect, as sources without a BB continuum are 
already fainter than others. We thus have three distinct 
groups in our sample. In Fig. |] we display the combined 
SWS spectra of the fourteen HAEBE objects arranged in 
these three groups. 

If we neglect the silicates, then the shape of the IR 
spectra of group la and lb sources is very alike. They both 
have a prominent cool dust component, rising in the mid- 
IR. What remains of the group lb spectra after continuum 
subtraction is essentially the same as for those of group 
la sources without their silicate features. This is shown 
in Fig. H where we compare HD179218 (group la) with 
HD100453 (group lb). Adding the right amount of amor- 
phous and crystalline silicate components, we can obtain 
the spectrum of HD179218 starting from HD100453. It is 
striking that both groups seem to have a similar contin- 
uum, yet the silicates behave completely different in both 
sources. 

To summarize these observational data, it first appears 
that both the near-IR and sub-mm excesses are similar for 
all stars in our sample; the mid-IR flux, on the contrary, 
shows large source to source differences. Furthermore, the 
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Table 2. Spectral slope of the power-law continuum fits 
in the IR region and the slope of the continuum in the 
sub-mm region. The slopes were measured in a log XF\ 
versus log A scale. For 51 Oph, we could not determine a 
power-law continuum in the IR. The temperature of the 
BBs needed to fit the remains of the power-law subtracted 
continuum are listed as well. HD 100453 and 51 Oph have 
no sub-mm slope listed because there are no 3 a sub-mm 
measurements available for these stars. For stars indicated 
with a '*', we had to include an IRAS flux (100 /mi flux 
for HD100546 and HD104237; 60 /zm flux for HD150193) 
in the determination of the sub-mm slope, as there is only 
one sub-mm measurement available for those stars. 



Group 


Source 


IR 


BB 


sub-mm 






slope 


Temp. (K) 


slope 


la 


AB Aur 


-1.20 


93 


-4.28 


la 


HD 100546 


-1.17 


170 


*-3.08 


la 


HD 142527 


-1.15 


73 


-3.60 


la 


HD179218 


-1.09 


195 


-4.28 


lb 


HD 100453 


-1.03 


148 




lb 


HD 135344 


-1.43 




-3.28 


lb 


HD139614 


-1.31 


169 


-3.20 


lb 


HD169142 


-1.25 




-3.24 


Ila 


HD104237 


-1.10 


absent 


*-2.81 


Ila 


HD 142666 


-1.03 


absent 


-2.91 


Ila 


HD 144432 


-1.32 


absent 


-2.65 


Ila 


HD150193 


-0.95 


absent 


*-2.67 


Ila 


HD163296 


-1.07 


absent 


-2.94 


Ila 


51 Oph 




absent 





spectra can be decomposed into at maximum three parts: 
a power-law, a BB and the solid state bands. 



spectral slope will be \F\ cx A~ 3 , while for small grains 
this slope will be equal to AF\ cx \~( 3+p \ where p is 
defined from Q u cx v p , it is the slope of the emissivity 
law of the dust at sub-mm wavelengths. Estimating the 
size based on spectral indices should be done with some 
caution, since the grain emissivity and the temperature 
distribution of the grains in the disc also affect this slope. 

Spectral sub-mm slopes have been determined for our 
sample and are shown in Table 0. Our sample shows a 
fair similarity in spectral slope, except for the stars AB 
Aur and HD179218. For AB Aur, the determination of the 
slope is more accurate than for HD179218, because the last 
source has only two sub-mm points at a small A-interval. 
As was already noted by van den Ancker et al. (1999), 
AB Aur has a very steep spectral slope [\F\ cx A -43 ); 
both stars jump out when compared to the average of our 
sample (cx A~ 32 ); the slope of AB Aur and HD179218 is 
significantly steeper than that of the tail of a black body 
(cx A -3 ). Bouwman et al. (2000) show that in AB Aur the 
10 /tm silicate band and the mid-IR continuum are domi- 
nated by micron-sized grains, while the millimetre contin- 
uum is produced by grains with typical sizes of several 100 
microns. These grains are optically thin at sub-millimetre 
wavelengths, contrary to millimetre-sized grains around 
other HAEBE stars (such as e.g. HD163296, van den An- 
cker et al. 2000). Possibly the grain size distribution in 
HD179218 is more similar to that of AB Aur, with some- 
what larger grains (but still in the range of ~ 100 /im size) 
producing the sub-millimetre flux. Sylvester et al. (1996) 
already noted from their mm/sub- mm photometry that 
the dust grains around most of their Vega-like systems are 
much larger than those found in the interstellar medium 
(sub-/tm sized, Mathis et al. 1977). 



2.2. Properties of the cold grains: mass and grain size 

The sub-millimetre flux is substantial in all our pro- 
gramme stars, with the exception perhaps of 51 Oph. This 
implies that a large amount of cold, large grains must be 
present further away from the star. The group-averaged 
sub-mm excesses are comparable. To compare the mass 
of the cold dust between different groups, we calculated 
Fgoo^mD 2 , a normalization of the cold dust mass which 
assumes that the stars are at the same distance. These 
data are also listed in Table |l|. The amount of cold dust 
does not differ substantially between the three groups. 

A first indication of the size of the cold grains can be 
obtained by inspecting the far-IR to sub-mm region. All 
our programme stars show a turnover in their SED at far- 
IR wavelengths, indicating that at sub-mm wavelengths 
we observe the Ray leigh- Jeans (R-J) SED of the ensemble 
of cold grains present in the disc. If these grains are large 
compared to the wavelength at which they radiate, the 



3. ISO-SWS spectra and their wealth of features 

We describe general trends in the appearance of the solid 
state bands we have detected in our sample of stars. We 
have carefully inspected the individual spectra in order to 
verify the reality of solid state bands. This can be done by 
analyzing the individual detector scans, and by inspect- 
ing the two independent scan directions at which the data 
were taken. We list the solid state bands and their possible 
identifications in Table [|. The features for which we have 
an unsure identification are listed with a question mark. In 
particular, the PAH feature at 11.2 /im can be easily con- 
fused with the olivine feature at 11.3 /tm. Caution should 
also be taken around 12 /tm, where an SWS band jump oc- 
curs. In Figs. ^ and [?], we show the continuum subtracted 
spectra in the 6 to 14 /tm wavelength region and the 15 
to 30 /im wavelength region respectively, highlighting the 
solid state bands found in these spectral regions. Below, 
we briefly discuss the different solid state components. 
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Table 3. Presence of the solid state and PAH bands. ■>/: detection, ?: possible detection, -: no detection. Regions for 
which we have no data or which we can not trust because of poor quality are indicated with n.a. (not available). 



Feature 


PAH 


PAH 


PAH 


PAH 


Si-O 


PAH" 


Oliv. 


Oliv. 


O-Si-O 


Oliv. 


FcO 


Oliv. 


Oliv. 


Oliv. 


H 2 ice 


A [/im] 


[3.3] 


[6.2] 


["7.7"] 


[8.6] 


[9.7] 


[11.2] 


[11.3] 


[16.3] 


[19] 


[19.8] 


[23] 


[23.8] 


[27.9] 


[33.7] 


[43.8] 


AB Aur 




V 


V 


V 


V 


V 






V 




V 










HD100546 


V 




V 


V 


V 


V 


V 


V 


V 




V 




V 


V 


V 


HD142527 


V 


V 






V 


7 


V 




V 




V 




7 


7 


v c 


HD179218 


V 


V 


V 


V 


V 


? 


v b 


7 


V 


V 


7 


V 


7 


%/ 




HD100453 


V 


V 


7 






V 




















HD135344 


? 




7 






7 












n.a. 


n.a. 


n.a. 


n.a. 


HD139614 
































HD169142 


V 


V 


V 


7 




V 




















HD104237 










V 




? 




V 


V 


? 


7 


? 


7 




HD142666 


V 


? 


V 


? 


V 


V 






V 






? 




7 




HD 144432 










V 




? 




V 




V 


n.a. 




n.a. 


n.a. 


HD150193 










V 




V 




V 




7 


n.a. 


n.a. 


n.a. 


n.a. 


HD163296 










V 




V 


V 


V 


V 


V 


? 


7 


7 


? 


51 Oph 










V 




? 




V 




7 


n.a. 


n.a. 


n.a. 


n.a. 



a possible blend with the olivine 11.3 [an band 
6 possible blend with the PAH 11.2 band 

c doublepeaked at 43.8 fim, probably also due to montmorillonite (Malfait et al. 1999a). 



3.1. The 2-7 micron wavelength region 

From Fig. [|, one can observe immediately that there is a 
striking similarity between the spectra from 2 to 7 ^m in 
all three groups, indicating a similar composition and 
temperature distribution of the material emitting in this 
range. This near-IR excess must be caused by small hot 
grains, most likely metallic Fe (or FeO), or a carbonaceous 
component, e.g. graphite or amorphous carbon (van den 
Ancker et al. 1999, Bouwman et al. 2000), as other mate- 
rials would not survive the high temperatures (<; 1000 K) 
close to the star. 



3.2. The PAH bands 

PAH bands are present in group la and lb sources, and 
only weakly in one source of group II (HD142666); they 
are seen in at least 50% of our total sample. This num- 
ber might still go up as weak features can be lost in the 
noise. PAHs are strongest in early type sources, as can be 
expected (since they are excited by UV radiation). That 
the PAHs most probably belong to the stellar environ- 
ment can be ascertained from e.g. ISOCAM data of HD 
179218, a source with is not extended (Siebenmorgen et 
al. 2000). Another interesting observation is their different 
appearance, at e.g. A ~ 6.2 /mi. PAHs in the ISM and in 
HII regions peak at A ~ 6.20 /im, while at our stars they 
rather peak at A ~ 6.25 /im. A preliminary result in this 
context is that non-extended sources tend to have PAHs 
at 6.25 /im (Van Kerckhoven, private communication). A 
more detailed study will follow (Van Kerckhoven C, in 
preparation; Peeters E., in preparation). 



3.3. The silicates 

In Fig. |^ and ^ we show the regions surrounding the 
strongest amorphous silicate bands. We also included a 
spectrum of the M supergiant /i Cep (Kemper, private 
communication) and the galactic center (Lutz et al. 1996, 
Tielens, private communication), as prototypes of amor- 
phous silicate. The spectra of group lb sources are not 
shown, as they do not show silicate bands. The 8-12 /tm sil- 
icate features are significantly different from the so-called 
'astronomical silicate', in the sense that they peak at 
11 /im rather than at 9.7 /tm. The only exception is AB 
Aur, a group la source. The peak shift is attributed to a 
change in the grain size distribution towards larger (/tm- 
sized) grains, and/or the presence of crystalline olivine, 
causing a peak at 11.3 /im (Bouwman et al., in prepara- 
tion) . 

In some sources such as HD100546, there is a large 
amount of crystalline silicates (Malfait et al. 1998a), while 
in other sources, such as HD163296, it is less so. That both 
AB Aur and HD100546 are members of group la, shows 
that the crystallization degree of the silicate material is 
independent of the shape of the overall spectrum. The in- 
terpretation from sub-mm data that AB Aur has the least 
processed dust (smallest grains at sub-mm wavelengths) is 
here further supported by the lack of crystalline features in 
its spectrum. We, however, derived a similar slope for the 
sub-mm region of HD 1792 18 (see Sect. 2.2), a source with 
a large amount of crystalline silicates. The determination 
of this slope was less accurate, but it is for sure steeper 
than the other sources. From this we can infer that co- 
agulation and crystallization processes occur on different 
time-scales. 
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Fig. 6. Continuum subtracted spectra of the 10 /im silicate feature. We also show the spectra of the M supergiant 
/i Cep and the galactic center for comparison. The galactic center spectrum has been converted from absorption to 
optical depth, which is plotted on the y-axis. Amorphous silicate, as seen in the Galactic Center spectrum, peaks 
at 9.7 /im, while crystalline silicates peak at longer wavelengths (~ 11 /im). Notice the quasi omnipresence of an 
11.3 /xm band, due to crystalline silicates and/or PAHs in the spectra of the Herbig Ae/Be stars. Only AB Aur has 
no crystalline signatures at all. Group lb sources are not shown due to a lack of band structure in this region. 



3.4- The 23 micron feature 

The [15-30] /im spectra (see Fig. 0) of group la and Ha 
sources are even more diverse than the [6-14] /xm spectra. 
From the spectrum of /i Cep, we can see that the amor- 
phous silicates alone can not account for the broad bands 
seen in other sources, e.g. AB Aur and HD142666. An ad- 
ditional component, around 22-28 /im must be present as 
well. This component can be attributed to FeO (van den 
Ancker et al. 2000) and/or crystalline silicates. 51 Oph 
shows two clearly separated bands, supporting our inter- 
pretation of the two component broad band. The same 
region for group lb sources can be fitted very well with 
black bodies with a temperature ranging between 150 and 
170 K. Also here we see no indication for silicate bands in 
group lb sources. 



3.5. Notes on individual sources 
Group la: 

AB Aur: The 10 /ion region of AB Aur shows a broad 
silicate feature superposed on a rising continuum. The 
feature peaks at 9.7 /im, which is typical for amorphous 
silicate. It is the only source in our sample to show merely 
pure amorphous silicates. The 20 /im region shows a 
broad band, probably due to a combination of silicates 
and iron oxide. PAH bands are moderately strong, and 
present at 6.2, 7.7, 8.6 and 11.2 /im. This star has the 
steepest spectral slope at sub-mm wavelengths. For a 
thorough analysis of this star, we refer to van den Ancker 
et al. (1999) and Bouwman et al. (2000). 
HD100546: This star is thoroughly analyzed by Malfait 
et al. 1998a. Amorphous silicate emission bands are 
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Fig. 7. Continuum subtracted spectra of the 15 to 30 /im region. We also show the spectra of the M supergiant (i Cep 
(amorphous silicates) and comet Hale-Bopp (crystalline silicates) for comparison. From a comparison with /i Cep, it 
is clear that most sources need an additional band longwards of the 18 /^msilicate feature to explain their spectra. 
Group lb sources are not shown due to a lack of band structure in this region. 



visible at 10 and 18 /an. It has the largest amount of 
crystallinity in our sample, and shows strong PAH bands 
at 3.3, 6.2, 7.7, 8.6 and 11.2 /im. The similarity with 
the SWS spectrum of comet Hale-Bopp is remarkable 
(Waelkens et al. 1998). The comet Hale-Bopp is an end 
result of grain processing in our solar system. Given the 
large similarities in spectral appearence with HD100546, 
we can assume that the same kind of dust processing has 
taken place in both sources. 

HD142527: The silicate feature at 10 /im is shifted to 
longer wavelengths, pointing to the presence of crystalline 
silicates. The continuum rises from longer wavelengths on 
than most sources, pointing to colder dust. Weak PAH 
bands are present at 3.3, 6.2 and 11.2 fim. We refer to 
Malfait et al. (1998a) for more details. This star is special 
because it has a large Lir/L* ratio, typical for stars 
with an active disc. But these stars also show evidence 
for outflow, what is not observed in HD142527. Another 



explanation for the large IR luminosity is the presence of 
a more embedded companion, but this still needs to be 
investigated. 

HD179218: The overall shape of this spectrum is quite 
similar to that of HD100546. It shows the 10 /im silicate 
feature, and has a rising continuum. The abundance 
of crystalline silicates is less, but still quite high. This 
source has a higher abundance of crystalline pyroxenes 
than crystalline olivines, unlike observed in other sources 
(Malfait 1999b). Strong PAH bands are present at 3.3, 
6.2, 7.7, 8.7 and 11.2 urn. 



Group lb: 

HD100453: This star does not show silicate emission 
bands. It has a rising continuum and moderately strong 
PAH bands are present at 3.3, 6.2, 7.7 and 11.2 /im. 
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HD135344: In this source the silicates are absent. 
The spectrum shows a rising continuum, but beyond 
27 fim the spectrum is no longer usable due to low S/N. 
Weak PAH bands are present at 3.3, 7.7 and 11.2 /im. 
HD 1396 14: Also here there are no silicate bands. The 
spectrum shows a rising featureless continuum, PAH 
bands are absent. 

HD169142: Moderately strong PAH bands are present 
at 3.3, 6.2, 7.7, 8.7 and 11.2 /im. The spectrum shows a 
rising spectrum without silicate emission bands. 



Group Ila: 

HD104237: The 10 /im spectrum shows a very strong 
silicate feature, superimposed on a flat continuum. The 
feature peaks longwards of 10 /im, so that crystalline sil- 
icates and/or larger silicate grains must be present. PAH 
bands are absent. 

HD142666: The spectrum shows a flat continuum, upon 
which the 10 /im feature is superimposed, peaking at 
10.3 /im. Very weak PAH bands are present at 3.3, 6.2, 
7.7, 8.6 and 11.2 /im. The 20 /im region shows a broad 
band, probably consisting of silicate and FeO. 
HD144432: The spectrum of HD144432 is very similar 
to that of HD142666, but PAH bands are absent. Inter- 
estingly, these 2 sources are in general quite similar, the 
most important difference being the inclination of their 
disc (Meeus et al. 1998). The 20 /im region shows a broad 
band, probably consisting of silicate and FeO. 
HD150193: This source has a very strong 10 /im feature, 
superimposed upon a flat continuum. The shape of the 
spectrum around 20 /im is very peculiar, and we are not 
sure if it is an artefact or real. PAH bands are absent. 
HD163296: The 10 /im silicate feature is superposed on 
a flat continuum. The 20 /im region shows a broad band, 
probably consisting of silicate and FeO. PAH bands might 
be present, but we cannot ascertain this. We refer to van 
den Ancker et al. (1999) and Bouwman et al. (2000) for 
more details and modelling. 

51 Oph: This source is the most extreme, in the sense that 
its turn-over point towards the R-J tail already starts at 
~ 30 /im; from which we can conclude is has a smaller 
amount of dust. It also has a very rich spectrum, with 
both gas (CO2 at 4.2 /im and H 2 at 5-6 /im) and solid 
state bands. The 10 /im silicate feature is strong and im- 
posed on a descending continuum. The 20-30 /im region 
shows two clearly separated bands, which we can attribute 
to amorphous silicate (18 /im), and to FeO or crystalline 
silicate (23 /im). The evolutionary status of 51 Oph is 
not very clear. Some authors classified this object as a Be 
star (Slettebak 1982), while others consider it as a Her- 
big Ae/Be star (Malfait et al. 1998b) or even a Vega-type 
star (Sylvester et al. 1996). A more detailed study of this 



object will be presented elsewhere (Meeus et al./van den 
Ancker et al., in preparation). 

4. Discussion 

Two important results of the comparative study of the 
SEDs of the programme stars concern the near-IR and 
the sub-mm excess: 1. a similar near-IR (1-8 /im) excess 
is observed for all stars. Small hot grains close to the 
star must be replenished because they are continuously 
destroyed by the UV radiation of the star. The overall 
IR spectrum being so diverse, but the near-IR so similar 
indicates that the material close to the star is homoge- 
nized; 2. the sub-mm excess is substantial for all stars in 
our sample. This implies that large grains already formed 
when the star formation process comes to an end, and that 
large grains remain present around the stars during their 
evolution towards Vega-type stars. In a survey of T Tauri 
Stars (TTS), Beckwith et al. (1990) also found that the 
disc mass does not decrease with increasing stellar age. 
Warm grains, however, seem to disappear on a shorter 
time-scale. 

The third remarkable observation is the strong varia- 
tion in strength of the silicate feature: for some stars it 
dominates the ISO spectrum, for others it is moderately 
strong, and for some stars it is even absent. It is surpris- 
ing that there is no relation between the silicate feature 
at 9.7 /im and the near-IR excess, although both emis- 
sion features must be caused by hot material, presumably 
located in the same region, and consisting of grains of 
a similar size. We can already exclude inclination angle 
effects since our sample includes two objects, HD142666 
and HD144432, with a very different inclination (Meeus 
et al. 1998) but with almost identical ISO spectra, both 
showing a prominent 10 /im silicate band. There are sev- 
eral possible explanations to our observations, and these 
'scenarios' will need to be confirmed by more detailed ob- 
servations and by careful modelling of the CS material. In 
the following subsection we propose a global model. 

4-1- Geometry of the disc and its effects on the SED 

As stated before, the main difference between group I 
and II sources is the amount of mid-IR excess, which is 
dominant and rising for group I sources, while moderate 
and rather descending for group II sources. Both groups 
also differ as far as the total IR luminosity is concerned: 
Lir/L» is on average 0.52 and 0.17 for group I and II re- 
spectively (see Table [l]) . Group II sources thus have the 
smallest emitting surface, and probably the smallest mass 
of warm dust. On the other hand, the solid state bands 
are present with equal average strength in group la and 
Ila sources: the solid state bands thus must be formed 
in yet another region. Natta et al. (2000) calculated the 
silicate 10 /im feature intensity with the models of Chi- 
ang & Goldreich (1997), and they needed to add a power- 
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law component to the emission of the disc atmosphere to 
fit their 10 micron spectra of T Tauri Stars; this finding 
supports our distinction between the region in which the 
solid state bands are formed and the region from which 
the power-law emission originates. 

The classification of our sample of HAEBE stars into 
two main groups could be explained with the following 
simple physical picture, which is shown schematically in 
Fig. H (upper panel: a group I source, lower panel: a group 
II source) and consists of the following components: I) a 
(partially) optically thin inner part (~ 10 AU in size); II) 
an optically thick, geometrically thin disc which forms in 
an early stage; and III) a flaring part, exposed to stellar 
radiation. Similar geometries were suggested in studies of 
protostellar accretion discs by e.g. Bell et al. (1997) and 
Nelson et al. (2000). There it is shown that the protostel- 
lar discs can have an inner disc (~ 10 AU in size) with 
increased scale height (i.e. "puffed up") which can shield 
the outer parts of the disc so that no flaring occurs. One 
can explain the SEDs of group II sources by assuming that 
the inner part (component I lower panel Fig. ||) is partially 
optically thick, shielding the outer parts of the disc from 
direct stellar radiation, which prevents the disc from flar- 
ing (hence one does not observe a BB component). SEDs 
of group I sources then can be explained by an additional 
flared component outwards of an optically thin inner disc 
which does not shield. The non-occurence of a group lib 
source could be a selection effect, since one would expect 
such sources to have the lowest flux levels. 

The emission of an optically thick, passive disc varies 
as AF\ oc A -13 (Friedjung 1985). The average IR slope for 
the sources in our sample is -1.2, very close to the value for 
the slope of optically thick discs. There is, however, some 
dispersion in our slopes, but nevertheless the assumption 
of an optically thick disc to explain the power-law com- 
ponent seems to hold when confronted with theoretical 
models. 

In what follows we discuss the separate components of 
our model and how they appear in the SED: 

— the geometrically thin, optically thick midplane 
consists of large grains, and is responsible for the IR 
power-law and the sub-mm continuum emission. Here 
resides the bulk of the material (Bouwman et al. 2000). 
The small hot grains in the inner part of the disc are 
responsible for the near-IR excess. The similarity of 
the spectra in this wavelength range suggests a very 
similar chemical composition and temperature distri- 
bution. 

— hot grains in a (partially) optically thin "disc atmo- 
sphere" cause the 10 and 18 /im amorphous silicate 
emission bands. That the solid state bands are indeed 
likely not to originate from the flared region around 
the disc is evidenced by: 1) our group Ha observations 
(solid state bands are present, but the BB caused by 
the flaring is absent); 2) the fact that silicate grains 



causing the 10 /im emission are too large to be tran- 
siently heated, so they have to be located close to the 
star; and 3) observations of evolved objects by Mol- 
ster et al. (2000), who notice that disc sources have a 
high abundance of crystalline silicates. Where this op- 
tically thin region is located exactly is not clear, but 
it certainly must be close to the disc and to the star. 
— the geometrically thick, optically thin, flared dust 
layer below and above the midplane is well mixed 
with the gas and contains small, warm grains that 
cause the 100-200 K black body component. Only a 
small amount of the total disc mass is located in this 
flared region. In group II sources this region must be 
very small or even absent, since we only see a power- 
law continuum in the IR. This conclusion is also sup- 
ported by theoretical models from Kenyon & Hart- 
mann (1987), which give an upper value for Lir/L* 
for flat discs of 0.25, while flaring discs have larger 
IR luminosities. When comparing the values we deter- 
mined for Lir/L» (see Table |l|) with their predictions, 
we see that group II sources indeed have Ltr/L* values 
arguing for a flat disc. The scale- height of the flaring 
determines the amount of mid-IR excess: the larger 
the flaring, the stronger the excess. The temperature 
of the BB is determined by the distance of the on- 
set of the flaring to the central star: the further away 
the flaring starts, the colder the BB. If we e.g. com- 
pare HD179218 with HD142527, then the flaring must 
start further away for the latter object, as evidenced 
by its BB temperature (see Table ||). That PAHs are 
only present in group I sources, with the exception of 
HD142666 (group Ha, but only very weak PAH fea- 
tures), suggests that PAHs are most likely located in 
the flared region, where they are exposed to the stellar 
UV radiation. They are small enough to be transiently 
heated, without needing to be located close to the star. 

Other geometries interpreting the SEDs of HAEBEs 
have been proposed, e.g. Miroshnichenko et al. (1999) pro- 
pose an envelope in addition to a disc to model the dust 
emission from HAEBEs. We note, however, that the highly 
abundant crystalline silicates in HD100546 (that have a 
temperature of 200 K or less, Malfait et al. 1998a) are 
most probably formed in the disc (Molster et al. 1999), 
suggesting that the grains at this temperature are not in 
a loosely bound envelope, but are intimately connected to 
the disc. We suggest that the 200 K black body compo- 
nent in this star, and by analogy in other HAEBE stars 
in our sample, is associated with the (flared) disc and not 
with an envelope. Waelkens et al. (1994) and van den An- 
cker et al. (1997) interpret the apparent broad dip around 
10 micron in the SEDs of HAEBEs as a physical gap in 
the radial distribution of the CS dust (there is simply no 
dust at a certain distance to the star, corresponding with 
the region where the 10 micron flux should come from). 
To cause a physical gap, another body surrounding the 
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central star must be present as well. In the light of these 
different possibilities, detailed spatial information is es- 
sential to disentangle the location of the different spectral 
components. 

4-2. Evidence for grain growth 

Evidence for grain growth in the discs of HAEBEs has 
been found by several authors. Radiative transfer mod- 
elling by e.g. Bouwman et al. (2000) shows that the dust 
grains around the Herbig Ae stars AB Aur and HD 163296 
are much larger than those of the ISM. Furtheron, Grady 
et al. (1996) observe accreting CS gas in HAEBEs and at- 
tribute this to large infalling objects. A population of large 
(~ 0.1-lmm) grains is needed to explain the observed sub- 
mm fluxes in our sample. These observations show that 
grains growth indeed takes place around HAEBE stars, 
and that it is an on-going process which we can observe 
indirectly by looking at a large sample of objects. 

We now consider the cause of the difference in SED 
between group I and II sources. If our interpretation of 
a flared region above and below the disc causing the BB 
component is correct, its absence in group II sources may 
imply that these small, warm grains have been removed 
due to coagulation and/or to radiation pressure exerted 
on the grains by the central star, so that the small warm 
grains have slowly disappeared and their scale-height has 
diminished accordingly. The absence of PAH bands in 
group II stars already supports the assumption that small 
grains are removed in the extended region. If we expect 
grains to grow during the star's evolution towards the MS, 
then the excess as a whole should decrease. This is con- 
sistent with the amount of IR luminosity we derive from 
the SEDs: stars from group II show a smaller Lir/L* ra- 
tio than group I stars, so that group II sources may have 
the most evolved dust grains. This simple evolutionary as- 
sumption is supported by observations of TTSs by Beck- 
with et al. (1990), where it is shown that older discs tend 
to be colder and less luminous. The grain-growth assump- 
tion to explain the differences between group I and II also 
holds in the sub-mm region: from Table || it is clear that 
group II sources have a less steep sub-mm slope than group 
I sources (on average -2.8 versus -3.6), which means that 
the latter have smaller grains radiating in this region. 

To conclude, the SEDs are consistent with a disc model 
in which the differences between group I and group II can 
be explained by a different extent of the warm flared layer 
above and below the mid-plane. Both the Lir/L* ratio 
and the sub-mm slope suggest that group II sources have 
larger grains than group I sources. These observations are 
consistent with an evolution from group I to group II 
sources. However, a larger sample and spatial information 
are needed to prove such an hypothesis. 



4-3. Influence of the stellar age on dust properties 

The age of the sample stars was derived by van den Ancker 
et al. (1998) using Hipparcos parallaxes, and is listed in 
Table [l]. Unfortunately, we do not dispose of ages for group 
lb stars. We are aware of the fact that our sample is biased 
towards the more evolved sources, since our sources are 
isolated. However, we can already reach some conclusions. 
There is no clear trend between age of the central star 
and amount of crystalline material: the star HD100546 
(showing the largest amount of crystalline dust) is proba- 
bly the most evolved one, while HD179218 (probably the 
youngest source) also shows a substantial amount of crys- 
talline dust. AB Aur, on the other hand, is also already 
more evolved, but shows only evidence for amorphous sil- 
icates. The range in age is very similar for group la and 
group Ha stars; from a confrontation between both groups 
we can conclude that the stellar age has no (or little) in- 
fluence upon: 1) coagulation, since group Ha stars have a 
less steep sub-mm slope than group la stars; 2) amount of 
material in the flared region, as this is much less or even 
absent in group Ha stars; and 3) presence of PAH bands, 
as they are merely absent in group Ha stars. We thus con- 
clude that the timescales on which star and disc evolve are 
not strongly coupled for the sample of HAEBEs studied 
here. 

4-4- The amorphous silicate behaviour 

It is surprising that the silicate feature is absent in group 
lb sources, as the rest of the SED argues for similar disc 
properties as for group la sources. It is not unreasonable 
to assume that grains grow during the stellar evolution 
towards the MS. Therefore, the absence of the silicate fea- 
ture could be easily explained by the absence of small 
grains (with average sizes less than a few /zm). However, 
the presence of a near-IR excess points to the presence of 
small, hot grains. Furthermore, almost all of the group lb 
sources show PAHs in their spectra, which are also caused 
by very small particles. It thus seems that the small sil- 
icate grains evolve differently than other small particles. 
The explanation for group lb sources could be either that 
the inner part does not exists (e.g. could be geometrically 
thin, optically thick), or that there are no small (^ 50 /im) 
silicate grains. 

Remarkably, in an atmospheric abundance analysis by 
the authors and another analysis by Dunkin et al. (1997), a 
silicon depletion around 3 out of 4 of the group lb HAEBE 
objects was revealed, while sources from group la and Ha 
were shown to have solar abundances. The photospheric 
silicon depletion for group lb stars may further support 
that around these stars, silicates behave differently. There 
are three possibilities to remove the small silicate grains 
selectively: 1) a composition effect; 2) a size-effect: the 
silicates are too large to be seen in that wavelength-region; 



G. Meeus et al.: ISO spectroscopy of 14 Herbig Ae/Be stars 



13 



or 3) aggregates with other materials. In what follows we 
will discuss these different effects. 

1) A composition effect: We expect the strength 
of the silicate emission to be related to (among others) 
the amount of silicates, as the emitting dust around our 
sources is optically thin. However, it is unrealistic to as- 
sume that the sources which do not show silicate emission 
have no silicates at all in their disc. This would mean 
that there were no silicates in the material from which 
the group lb stars are formed. This hypothesis is most 
unlikely, since the material from which stars are formed 
is the ISM, which is relatively uniform in composition. 
There is no reason to assume that the silicates were not 
there in the beginning around some stars and were there 
around others. Besides, the presence or absence of silicate 
emission does not depend on the initial composition: two 
very young HAEBE stars, LkH Q 224 and LkH Q 225 are 
located close to one another, so must be formed out of the 
same material. Surprisingly, the first object does not show 
any silicate band, while the second object does show sili- 
cate absorption (van den Ancker 1999). This observation 
favours a scenario where the presence or absence of silicate 
emission is determined by other characteristics than only 
the chemical composition of the CS disc. 

2) A size effect: Since we found evidence for hot grains 
(from the 2-10 /xm excess) capable of producing a strong 
10 /xm feature if they are (partially) composed of silicates, 
we must come to the conclusion that the silicate grains 
around e.g. HD100453 and HD169142 (both group lb ob- 
jects) do not produce a silicate bump because these grains 
are on average larger than the wavelength (i.e. 10 /xm). 
Large silicate grains result in an emission resembling a 
black body without strong spectral signatures reflecting 
the chemical composition. Also Hanner et al. (1994) con- 
clude that the absence of small silicate grains is the cause 
for the weak silicate emission features in comets Austin 
and Okazaki-Levy-Rudenko. 

The absence of the 10 xtm silicate feature does not nec- 
essarily mean that also the 18 /jm silicate feature must be 
absent: e.g. in NGC 6302 (Lim ct al. 2000), the 10 /xm sil- 
icate feature is supressed because of the cold tempera- 
ture of the silicate dust and the dominant emission of 
the C-rich dust, but a 18 /xm feature is observed. There- 
fore, we searched the group lb sources for silicate fea- 
tures at 18 /xm. But unlike group la sources, group lb 
sources do not require solid state bands in addition to the 
BB+power-law to fit their continuum. We cannot fully ex- 
clude the possibility that small silicates are present in the 
discs of group lb stars, but they must be so by a much 
smaller amount and/or colder than in group la and Ha 
objects. Sylvester et al. (1996), however, claim to detect 
silicates (around 18 /xm) in the spectra of HD169142 and 
HD135344, indicating that these stars do have larger sili- 
cate particles; from our data, we cannot confirm this ob- 
servation, however. Modelling should determine how much 



silicate material can be present in the dust without being 
revealed in the spectrum. 

3) Another possibility is that, after coagulation, 
small silicate particles are locked up into larger grains, 
composed of both small silicate grains and some other ma- 
terial (as is seen in interplanetary dust particles (IDPs)). 
This would make them invisible if the mantle surrounding 
the silicate is sufficiently thick. If these coagulated grains 
are transported towards the star, the silicate material will 
start to evaporate when close enough, while Fe or carbona- 
ceous material can survive at higher temperatures. This 
scenario can account for both the presence of a near-IR 
excess and for the absence of hot silicate grains. 

5. Conclusion 

ISO-SWS spectra have shown that there is a large diver- 
sity concerning IR spectral features and shapes in Her- 
big Ae/Be stars. The results from this ISO-sample can be 
summarized as follows: 

1. our sample of 14 HAEBE stars can be classified into 
two main groups, based upon the shape of the contin- 
uum (flat or rising); a division which is further sup- 
ported by the decomposition of the continuum into a 
power-law and a BB 

2. the disc geometry is as follows: a geometrically thin 
disc being responsible for the power-law continuum, 
and a flared region of warm dust around the thin disc, 
causing a rising continuum in the IR. An optically thin 
inner disc causes the solid state bands. Hot dust in the 
inner part of the disc is responsible for the near-IR 
excess. This disc geometry has also been proposed for 
TTS (Chiang & Goldreich, 1997) 

3. the near-IR spectral region is very similar and the sub- 
mm emission substantial around all sample sources, 
indicating homogeneity of the hot material and sur- 
vival of large grains throughout the whole prc-main 
sequence evolution towards the MS 

4. the mid to far-IR region, on the contrary, is very di- 
verse, and we attribute this to the amount of flaring 
in the disc 

5. group I sources may evolve into group II sources; the 
latter have evidence for larger grains and lack the 
flared region present in group I sources 

6. the presence of PAH bands cannot be correlated to 
any of the stellar parameters, but they are only present 
in stars with a large amount of warm dust (group I). 
They are most probably located in an extended region 
around the disc, where they are irradiated by the star. 
Our PAH bands differ from PAH bands in the ISM 

7. surprising is the independent behaviour of the sil- 
icate grains: although other small particles are still 
present, small silicate grains seem to be absent around 
several stars. This poses an intriguing problem: what 
happens to the silicates, what causes them to behave 
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so differently? It is important to get a definitive an- 
swer on the warm silicates, to know to which extent 
they can be hidden. Only then can the observations 
and proposed model converge to a consistent picture. 
Therefore, the next step in our study will be a detailed 
modelling of some of the sources 
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Fig. 4. ISO spectra of the fourteen sample HAEBE 
stars, ordered by group. Group la: AB Aur, HD100546, 
HD142527, HD179218; group lb: HD100453, HD135344, 
HD139614 and HD169142; and group Ha: HD104237, 
HD142666, HD144432, HD150193, HD163296 and 51 Oph. 




Fig. 8. Schematic presentation of the model. The disc con- 
sists out of three parts: I) a (partially) optically thin inner 
part (~ 10 AU); II) a geometrically thin, optically thick 
midplanc; and III) a flaring part. The upper panel shows 
an entirely optically thin inner disc, where the star il- 
luminates the disc's surface and causes the flaring. The 
lower panel shows the case where the surface of the disc is 
shielded from direct stellar radiation by an optically thick 
region, so that it does not flare (Bell et al. 1997, Nelson 
et al. 2000). 



